The Statistics of the Work Done on a Quantum Critical System by Quenching a 

Control Parameter. 



Alessandro Silva 
Abdus Salam ICTP, Strada Costiera 11, 34100 Trieste, Italy 
(Dated: September 29, 2008) 

We study the statistics of the work done on a quantum critical system by quenching a control 
parameter in the Hamiltonian. We elucidate the relation between the probability distribution of the 
work and the Loschmidt echo, a quantity emerging usually in the context of dephasing. Using this 
connection we characterize the statistics of the work done on a quantum Ising chain by quenching 
locally or globally the transverse field. We show that for local quenches starting at criticality the 
probability distribution of the work displays an interesting edge singularity. 
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A series of recent experiments with cold atomic 
gases [H Q spurred new interest on the dynamics of quan- 
tum correlated systems. A number of fundamental issues 
on the nonequilibrium physics of quantum systems are 
being addressed, ranging from the relation between non- 
integrability and thermalization Q , to the universality of 
defect production for adiabatic quenches across quantum 
critical points jij. In this broad context, a paradigmatic 
example of experimental protocol is the instantaneous 
quench: an abrupt change, either global or local, of a 
control parameter g from some initial value g to a final 
one g±. Experimentally, it has been shown that the dy- 
namics after such quenches may show intriguing features, 
such as collapse and revivals of the order parameter for 
quenches done through a quantum critical point [j], [B| , as 
well as the absence of thermalization in systems close to 
integrability 0, 

Theoretically, the study of quantum quenches received 
considerable interest: after a series of classic works on the 
nonequilibrium dynamics of the quantum Ising model Q , 
recent investigations focused on characterizing the long 
time asymptotics of correlation functions 0, [1] , their be- 
havior as compared to their thermal counterparts [|[ , and 
the universality emerging in the quench dynamics at a 
quantum critical point [9(. Partial information on the 
internal dynamics of the system can be obtained in a 
variety of ways. One may extract the way in which exci- 
tations propagate by looking at the time dependence of 
correlators after a quench [a, [§| . More subtle informa- 
tion on the establishment of quantum correlations can 
be obtained by studying the dynamics of entanglement 
entropies [l(J. The purpose of this Letter is to discuss 
a more basic way to characterize both the internal dy- 
namics and the quench protocol itself by obtaining infor- 
mation on how far from equilibrium the system has been 
taken. This can be done by studying the statistics of a 
fundamental quantity: the work W done on the system 
by changing its parameters. 

The main observation behind this proposal is that the 
quench protocol resembles a standard thermodynamic 
transformation. However, since a quench takes the sys- 



tem out of equilibrium, the work W, unlike in a qua- 
sistatic process, is characterized by a probability distri- 
bution P(W) [HE El- Below, we focus on the char- 
acteristic function of P(W), defined as 



G(t) 



dwe wt p{w), 



(1) 



and study it for the prototypical example of quantum 
critical system, the quantum Ising chain. We first elu- 
cidate a useful relation between Git) and the Loschmidt 
echo, a quantity emerging in various physical contexts, 
most notably the Fermi edge singularity |l 4ll. q uantum 
chaos [15[ , and the physics of dephasing [lfj, ll7| . Using 
this connection and a combination of field theoretic tools, 
we compute exactly and analytically G(t) for global and 
local quantum quenches. In both cases, we characterize 
the fluctuations of the work and its probability distri- 
bution. Interestingly, we show that for a local quench 
starting at the quantum critical point the function P(W) 
displays an edge singularity. 

Let us start by briefly discussing the relation between 
the characteristic function G(t) and the Loschmidt echo. 
For a generic quench H(go) — > H(gi), the Loschmidt 
echo is defined as C(t) —\ Q(t) | 2 , where the amplitude Q 
is given by 



Q{t) 
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Here H(go) and H(gi) are the initial and final Hamil- 
tonian respectively, and the average is taken with re- 
spect to the initial equilibrium density matrix po = 
exp[— pH(go))/Z. The Loschmidt echo can be seen as 
a measure of the sensitivity of the system to the quench. 
The connection with P(W) emerges by noticing that for 
a generic quench the characterization of the work done on 
the system requires two energy measurements: one before 
and one after it 12, 13 1 . If the results of such measure- 
ments are E and E, the work done is then W = E — E. 
Hence if | <fr n ) are the eigenstates of energy E n of H(gi), 
we have that 



p ( w ) = 5 ( w - ( En - E m )) i <$„ 



Pm, (3) 
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where | are the eigenstates of Hq with energy E m , 
and P m = exp[— [3E m ]/Z. The characteristic function is 
then G(t) = ^ n e*(«»-Bm)* | ($„ | * m ) | 2 P m , which 
is readily recognized to be the complex conjugate of the 
amplitude defining the Loschmidt echo G(t) — [Q(t)]*. 
This equality is actually a special case of the generalized 
quantum Jarzynski equality [llL [l2| recently derived in 
Ref. llS for problems in which g is taken from go to g\ 
in a finite time interval along a generic path g(t). The 
Loschmidt echo C(t) can be in principle measured by 
studying the dephasing of an auxiliary two level system 
coupled to the system of interest [2, [la, [T3] ■ In the same 
setup, the probability distribution P(W) can be directly 
extracted from the absorbtion spectra associated to op- 
tical transitions in the auxiliary two level system (l8| . 

For a global quench the work done is extensive. There- 
fore in the thermodynamic limit the probability distribu- 
tion P{w) of the work per unit volume w = W/V will 
be a strongly peaked function, with fluctuations scaling 
as l/W. This suggests that P is a nontrivial function 
only for small systems or for local quenches. Despite this 
fact, it is interesting and instructive to study the work 
statistics for a paradigmatic example: a zero tempera- 
ture global quench of the transverse field in a quantum 
Ising chain [19j . The latter is defined by the Hamiltonian 



(4) 



where tr^' z are the spin operators at site i, J is an overall 
energy scale (below we set J = 1) and g is the strength of 
the transverse field. The one dimensional quantum Ising 
model is the prototypical, exactly solvable example of a 
quantum phase transition, with a quantum critical point 
at g c — 1 separating two mutually dual gapped phases, a 
quantum paramagnetic one (g > g c ) and a ferromagnetic 
one (g < g c ). 

Let us now consider a global change at time t = Q of the 
transverse field from an initial value go to a final one g± . 
The analysis of the Loschmidt echo can be efficiently per- 
formed after a Jordan- Wigner transformation [3] . In the 
fermionic representation, the Hamiltonian Eq.(j4]) takes 
the simple form 

k>0 

+ ism(k){c' k c[_ k - c-kCk), (5) 



where Ck are fermionic operators. The diagonal form 
H = Y,k>o E k(lhk ~ 7-fc7l fc )> with energies E k = 
^/(g — cos(fc)) 2 + sin(fc) 2 , is achieved after a Bogoliubov 

rotation c k = u k (g)jk - Wfc(ff)7_fc> c'_ k = u k {g)T_ k - 
iv k (g)"f k . The coefficients are given by 

u k (g) = cos(0 fc (») v k (g) = sm(9 k (g)), (6) 

where ta.n(28 k (g)) = sm(k)/(g — cos(fc)). In this rep- 
resentation, the Loschmidt echo for the quantum Ising 



model following both a global and a local quench of g has 
been recently shown to the expressible in terms of ma- 
trix determinants, which were afterwards analyzed nu- 
merically [3, [13] ■ Below, we compute analytically the 
Loschmidt echo employing field theoretic tools, which, in 
contrast to previous approaches, have the advantage of 
giving a clear insight on the physics of the problem. 

Our first task is to express the ground state | ^o) 01 
energy Eg of the initial Hamiltonian H{go) in terms of 
the eigenmodes 7fe diagonalizing H(g{). If we call rj k the 
eigenmodes of Ho it is easy to see that r\ k — Uk~1k ~ 
*Vfe7l fe , where 



Uk = u k (g )u k (gi) - 
V k = u k (g )vk(gi) 



- Vk(go)uk(gi)- 



(7) 

(8) 



Hence the equation r] k \ ^o) — 0, characterizing our ini- 
tial state can be easily solved giving 



1 



exp 
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v k t t 
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where Af is the normalization constant, and | 0) is the 
vacuum of the fermions jk- The structure of this state 
closely resembles that of integrable boundary states en- 
countered in statistical field theory. In particular, the 
amplitude Q is given by 



Q(t) =e lEot (^ I e' lH ^)t | *„) 



(10) 



-iSEt 
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where B{k) — —iVk/Uk, and SE — E\ — Eq, where Ei 
is the ground state energy of H(gi). Up to an irrelevant 
prefactor, Eq. lfTOf maps after a Wick rotation it — > r 
onto the partition function of a two dimensional classical 
Ising model constrained on a cylinder of height r and 
with boundary conditions on the two ends described by 
| ^o)- Hence, using techniques originally developed for 
integrable boundary states [201 ] it is easy to obtain 



G(t) 



-iSEt e L Jo 



■ log 
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(11) 



where L is the linear size of the chain. 

We can now derive all the cumulants C n of the proba- 
bility distribution of the work per unit length P(w) by ex- 
panding in power series log(G(t/L)) = J2n=o C n (it) n jn\. 
As expected, from Eq. lfTTj) we obtain that as L grows, 
C n oc l/L™ -1 , leading to the suppression of all fluctua- 
tions in the thermodynamic limit. In order to study the 
effects associated to the presence of the quantum critical 
point, let us focus on the dependence on <?o and g\ of the 
average work per unit length and of its variance. The 
first is in general given by (w) = SE / L + w e , where the 
excess work is w e > 0, in agreement with standard ther- 
modynamic relations. The latter and the variance have 



3 



0.5 
0.4 
0.3 
0.2 
. 1 



0.5 



1.5 
do 



2 . 5 



FIG. 1: A plot of the average excess work per unit length 
w e /gi (full line) and of the variance L/gl((Aw) 2 ) (dashed 
line) vs. go for a global quench from go to g\ S> 1. The 
presence of the quantum critical point at go — 1 is signaled 
by a discontinuity in the derivative. 



a particularly simple form if the final transverse field is 
large g i » I , in which case w e =3i(l+3o — I .9o~ 1 \)/^9o 
and((A U ;) 2 ) = (g 2 )/L(g$+4g i -3- S ign[g 2 -l}(g*-4g 2 + 
3))/8<7q. These functions are plotted in Fig.[TJ where one 
may see that both w e and ((Aw) 2 } signal the presence of 
the quantum critical point with discontinuities in their 
derivative at g = 1. More striking universal effects as- 
sociated to the quantum critical point are observed by 
studying the asymptotics of Q(t) at long times, governed 
by the long wavelength modes. If for example one looks 
at 9o^ 9i 7^ 1 one may easily obtain 



g(t) ~ g°° ex P 



A m\L 



( fn\ 



m 



2 g — 1irfi\t 

(imt)i 



, (12) 



where Q°° is the asymptotic value attained by Q (sig- 
naling a delta function peak in P(w)), A is a constant 
and uii =| gi — 1 |. Passing again to imaginary time, 
the dependence on t in the exponential corresponds to 
the dependence on the thickness of the free energy of the 
Ising model on a cylinder, which away from criticality is 
exponentially cutoff by the correlation length £ = I jm\ . 
At criticality, of course, it becomes a power-law as a re- 
sult of the establishment of long range correlations. A 
more detailed study of the statistics of P(W) for global 
quenches will be reported elsewhere 

Let us now pass to a much more interesting situation 
in which we expect the work done on the system to show 
nontrivial fluctuations: a local quench of the Hamiltonian 
from Hq — H(g) to H + V, where 



V = -6ga z . 



(13) 



In order to capture the main differences with the previous 
case, let us start by considering the case Sg <C 1 and 
evaluate Q(t) within a second order cumulant expansion 

0( f ) = (jHot e -i(H +V)t\ _ (r e -if* dt> V(t')^ 

= e -Hv)t e ~Uo dt i dt 2( T lv(ti)v(t 2 )])^ , 1A s 



where V(t) — exp[iH t]V exp[— iH t}. Using the 
fermionic representation of the spin operators we obtain 

S_9 
L 



V 



k,k' 



c k c\, 



c\c k , 



(15) 



Hence writing V in terms of the eigenmodes r\ k of Ho and 
substituting in Eq. (fT4"|) . with the help of Wick theorem 
we obtain 



e -iSBt e -f(t) 



(16) 



Here the energy shift SE is given by 



SE 



Sg 



dk 
2^ 



cos(29 k (g)) 



(Sg) 2 r dkdk' V(k,k') 



2 J_ n (2ir) 2 E k + E' k 



where V(k,k') = sm(29 k ) sm(29 k >) + 4cos(6» fe ) 2 sin(^) 2 . 

The most important information on the statistics of 
the work done on the system is contained in 



/(*) = ■ 



(Sg) 2 r dkdk' V(k,k r ) 



„(2n) 2 (E k + E' k 



(i _ e -i(E k +E k ,) t y 17) 



From this expression we may again estimate the various 
cumulants of P(W) by expanding in power series f(t). 
In particular, the variance is given close to the critical 
point by 



(aiT r)=(g)(2(1 ,r) I.H2-! lo, 



1 



8 



)) 



Despite the fact that this function has a logarithmic sin- 
gularity of the first derivative at g = I, as originally found 
in studies of dephasing 17| , it is important to notice that 
the integral leading to this expression gets contributions 
from all frequencies (not just small k). Hence universality 
does not emerge substantially. 

In order to obtain information on universal ef- 
fects one has to study the asymptotics of Q(t) for 
long times. This can be done by looking at the 
asymptotic value attained by / at infinity = 
(Sg) 2 /2fdkdk'/(2Tr) 2 V(k 1 k')/(E k + E k ,) 2 . Close to the 
quantum critical point g — 1, this is given by 



(18) 



'Sg' 


2 

log 


f 


_2tt_ 




Us-i|J 



Hence as t — > +oo we have 



£(t)~\g-l | 2 (^) 2 



(19) 



The Loschmidt echo vanishes at the quantum critical 
point with a cusp singularity. As shown below, the van- 
ishing of the Loschmidt echo is the result of an orthog- 
onality catastrophe, originating from the changing of a 
local scattering potential in a nontrivial, yet gapless, ef- 
fective fermionic system. In particular, if we set g = 1 
the long time decay of Q is a power law 



G(t) 



-iSEt 



(it)-(*) 



(20) 
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Hence we expect the probability distribution P(W) to 
display an edge singularity 



P(W) sa 6(W - 8E)(W - 5E)(%) 



(21) 

This expectation is readily confirmed by the exact so- 
lution of the problem for a local quench starting at the 
critical point g = 1. This can be obtained employing the 
scaling limit of the quantum Ising model in the Majorana 
representation 



Hq [ip, ip] = / dr i mip ip 



-If d r f + -f d r If, (22) 



where if and (p are Majorana fermions, and the mass is 
related to the transverse field by m = g — 1. The quench 
consists in going from Ho to Ha + V, where V[ip, if] = 
i 5m tp(0) <p(0). 

In order to compute Q at criticality (m = 0) let us use 
a trick due to Itzykson and Zuber [21| . We start by com- 
puting [(?(i)] 2 . Introducing two copies of the Majorana 
fermions, f\ t 2 and fi,2, we have 



[Q{t)f 



e Mot e 



-i(H +V)t 



(23) 



where Hq = H [ipi, ip~i] +H [ip 2 , $2] and V = V[fi,ip{\ + 
V[f>2, tp~2\- The most elegant way to proceed consists now 
in combining the Majorana fermions into Dirac fermions 
= (<Pi +i<P2)/V2 and ^> l = (<pi +^2) A/2, and then 
introducing a pair of nonlocal operators [22j defined as 
*+(r) = (*«(r) + * L (-r))/V2, and *_(r) - (* B (r) - 
^ r i,(— r))/\/2i. In terms of these 

Ho = J dr^ + {-id r )^+ + ^_{-idr)^-, (24) 

V = Jm(^(G)*+(0) - #t(0)*_(0)). (25) 

Physically, it is now evident that we have two chiral 
modes subject to local potential scattering of opposite 
sign characterized by phase shifts ±5 — ±5m/2. The 
computation of Q 2 is now a standard problem solvable 
by bosonization [3, [23| . We find that 



5(t) 



1 



1 + it 



(26) 



The complex conjugate of this expression is readily rec- 
ognized to be the characteristic function of the Gamma 
probability distribution 



P(w) 



w 



r 



ay 



(27) 



which indeed displays an edge singularity with a expo- 
nent (S/tt) 2 = (5g/(2Tr)) 2 consistent with the one ob- 
tained by the cumulant expansion. 



In conclusion, after elucidating the connection between 
the probability distribution P(W) of the work done on 
a system in a quantum quench and the Loschmidt echo, 
we characterized P(W) for global and local quenches of 
the transverse field in a quantum Ising chain. As men- 
tioned before, the experimental measurement of P(W) 
requires the realization of an optical absorbtion exper- 
iment in a fully controllable setting. Recent proposals 
for the realization of quantum spin chains using bosonic 
atoms in optical lattices [3] give a possible, concrete way 
to pursue this goal with the available experimental tools. 

I would like to thank G. Mussardo for important 
discussions, help and encouragement throughout this 
project and R. Fazio for stimulating my interest on the 
Loschmidt echo. I would also like to thank V. Kravtsov, 
A. Nersesyan, D.Rossini and G. Santoro for discussions. 



[1 
[2 
[3 

[5 
[6 

[-: 

[9 

[10 

[11 
[12 

[13. 
[14; 
[IS 

[16 
[17 
[18 
[19 

[20 

[21 

[22 

[23 
[24 



M. Greiner et al, Nature 415, 39 (2002); M. Greiner et 
ah, Nature 419, 51 (2002). 

T. Kinoshita, T. Wenger, and D. S. Weiss, Nature 440, 
900 (2006). 

M. Rigol et al, Phys. Rev. Lett. 98, 050405 (2007) 
W. H. Zurek, U. Dorner, P. Zoller, Phys. Rev. Lett. 
95 105701 (2005); A. Polkovnikov, Phys. Rev. B 72, 
161201(R) (2005). 

E. Altman and A. Auerbach, Phys. Rev. Lett. 89, 250404 
(2002). 

P. Mazur, Physica 43, 533 (1969); E. Barouch, B. M. 
McCoy and M. Dresden, Phys. Rev. A 2, 1075 (1970). 
K. Sengupta, S. Powell, and S. Sachdev, Phys. Rev. A 
69, 053616 (2004). 

F. Igloi and H. Rieger, Phys. Rev. Lett 85, 3233 (2000). 
P. Calabrese and J. Cardy, Phys. Rev. Lett, 96, 136801 
(2006). 

P. Calabrese and J. Cardy, J. Stat. Mech. (2007) P10004. 

C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997). 
J. Kurchan, e-print |arXiv:cond-m at/0007360v2. 

P. Talkner, E. Lutz, and P. Hanggi, Phys. Rev. E 75, 
050102 (2007). 

K. D. Schotte and U. Schotte, Phys. Rev. 182, 479 
(1969). 

R. A. Jalabert and H. M. Pastawski, Phys. Rev. Lett. 
86, 2490 (2001); Z. P. Karkuszewski, C. Jarzynski, and 
W. Zurek, Phys. Rev. Lett. 89, 170405 (2002). 
H. T. Quan, et al. Phys. Rev. Lett. 96, 140604 (2006). 

D. Rossini, et al. Phys. Rev. A 75, 032333 (2007). 
A. Silva, to be published. 

S, Sachdev, Quantum Phase Transitions , (Cambridge 
University Press, Cambridge, 1999). 
A. LeClair, G. Mussardo, H. Saleur, S. Skorik, Nucl.Phys. 
B453 (1995) 581-618. 

J. B. Zuber and C. Itzykson, Phys. Rev. D 15, 2875 
(1977). 

C. L. Kane, K. A. Matveev, and L. I. Glazman, Phys. 
Rev. B 49, 2253 (1994). 

Here we set the ultraviolet cutoff a ~ J = 1. 

L. M. Duan, E. Demler, and M. D. Lukin, Phys. Rev. 

Lett. 91, 090402 (2003). 



